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__i . Positron binding to anions is investigated within the work formalism proposed by Harbola and 

^ , Sahni (HS) for the halide anions and the systems Li~ through 0~ excluding Be~ and A'^". The 

, total ground state energies of the anion-positron bound systems are empirically found to be an upper 

fvg ^ bound to the Hartree-Fock(HF) energies. The computed expectation values as well as positron 

and positronium afBnities are in good agreement with their restricted Hartree-Fock counterparts. 

i-Q ' Binding of a positron to neutral atomic species is also investigated using an iterative method. 
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^ . I. INTRODUCTION 

CJ 

^^, Despite extensive studies [1-11] the problem of a positron forming a transient bound state with a neutral atom 



Oh, 



has not yet been decisively settled. It is known that a negatively charged ion always binds a positron; in fact it 
turns out that there is an infinite number of bound states. The classical Coulomb interaction between an anion 
and a positron is sufficient to bind a positron while in case of neutral atoms the polarization of electron charge 
distribution is found to be vital for binding |l2]. Schrader et al. [p3|jl^] have calculated the positronium affinities for 



> 

^^ I halogen atoms using a diffusion quantum Monte Carlo technique in which the core region is represented by a model 

vQ ', potential. Accurate calculations on the positronium hydride, positronium helium and positronium lithium systems 

^^ using the all particle diffusion Monte Carlo method have been recently reported by Yoshida et al. |Q,y . Very recently, 

^\ Bressanini ||l^ et al. have investigated the positron-anion bound states for the anionic species Li~ ,B^ , C^ , 0~ and 

c/3 , F^ using variational and diffusion quantum Monte Carlo techniques. They found that except B~ all these systems are 

" ^ I stable against dissociation into the corresponding neutral atom and positronium. Such methods are computationally 

^~^| very demanding thereby restricting their applicability to systems with fewer particles. The independent particle 

J^ ... 
O approximation methods such as the Hartree-Fock or density functional theories enable one to compute with case the 

^ ] anionic-positron interactions for all the anions in the periodic table. Calculations within the restricted HP theory 
k>( ' have been reported by Cade and Farazdel ||l6|,0 and Patrick and Cade |1^ for the systems Li~ through O^ as well 

JH , as for halogen anions. They have also examined the stability of anion-positron bound systems with respect to their 
dissociation into the corresponding atom and positronium. For the description of many electron-positron bound system 
within the density functional theory (DFT), Chakraborty [1^ has proposed a two component density functional theory 
wherein the electron and the positron densities play the role of basic variables. This was subsequently extended by 



Harrison |19| to incorporate self- interaction correction |gO[ wherein the calculated positron and positronium affinities 
are in close agreement with their HF counterparts. By invoking the Slater transition state concept Baruah et al. pl| ] 
have obtained decent estimates of positron affinities to negative ions. Earlier, density functional calculations of the 
positron-anion bound state had also been carried out by Kanhere et al. 123] . 

Recently, an attractive alternative to the Hartree-Fock theory was proposed by Harbola and Sahni p3. The 
exchange-correlation potential in the so-termed work formalism of Harbola-Sahni (HS) is obtained as the work done 
in bringing an electron in the electric field of its Fermi-Coulomb hole charge distribution. The work formalism HF 

1 



approximation essentially emerges if one ignores Coulomb correlations and the correlation contribution to the kinetic 
energy in the Kohn-Sham theory |24|. Recently, Holas and March ||2^ have shown that the HS potential can be 
obtained from the full second order density matrix. The total atomic energies and various one electron properties 
of work formalism are practically equivalent to those of HF theory [g6|. Unlike other local exchange-only density 
functionals, the work formalism of HS selectively gives convergent orbitals and eigenvalues for negative ions that are 



comparable to the HF accuracy 1 30 1 . In the present communication, we address the problem of positron binding to 
anions and atoms within the exchange-only work formalism. The purpose of the present work is two- fold: 1) To 
compute the positron affinities and binding energies using local orbital-independent density functional theory. 2) To 
test the work formalism of the Harbola-Sahni for the description of many electron-positron system. The computed 
positron affinities, binding energies and < ?'" > — moments (n = — 1 through 2) within the work formalism are 
compared against their restricted Hartree-Fock (RHF) counterparts. For the neutral atom-positron bound states we 
use an iterative method similar to the one used by Patrick and Cade [Q. The plan on the presentation is as follows: 
In section II we outline the theory of work formalism for the description of electron-positron system while section III 
deals with the results and discussion which will be followed by conclusions in section IV. 

II. THEORY 

The non-relativistic Hamiltonian H for an N-electron and one positron system is the sum of the electronic part H^ 
(atomic units are used throughout), 

consisting of the kinetic energy of electrons, the electron-nuclear interaction, the electron-electron repulsion, and the 
Hamiltonian Hp for a positron, 

H --^ ^ V ^ 

containing the positron kinetic energy, the nucleus-positron repulsion and the positron-electron attractive interaction. 
Within the exchange-only work formalism of Harbola-Sahni the local exchange potential is obtained as the work 
done in moving an electron in the electric field £ of its Fermi hole px(r,f') charge distribution, that is. 



W^{r) = - I £.,{r')-dl' (3) 

with 

^-(^1= / r7'.,:, i^~nd'r' (4) 

and 



Px{r,r') 1^ ^3 , 



'^^'^" y:;W)^^ ■ ^^ 

Now, in order that the effective potential experienced by electron is well defined the curl of the "exchange electric 
field" represented by Eq. (0) should vanish. This is the case for the closed shell atoms and open shell atoms in the 
central field approximation in which the present calculations have been carried out (see references [E3, p6[for details). 
It is to be the noted that HS exchange potential obtained this way differ from the exact Kohn-Sham potential only 
by the kinetic correlation contribution [E^pS. The orbitals, 0i, in Eq. (pi) are the solutions of the (Kohn-Sham like) 
HS equation 



Here, n^if) = ^,- | ipiif) p is the electron density, n+(r) =| ■0+(r) p is the positron density with i/j^ being the 
solution of the corresponding differential equation for the positron 

{-^ + f - / f^^d^^'Y^if) = e^i^^ir). (7) 

The effective potential seen by the positron in the exchange-only forinalism also has the interpretation as the work 
done in moving the positron in the field of the electronic and nuclear charge distribution. Eqs. (|6|) and (|7|) are solved 
self consistently to obtain the ground state energy of the electron-positron combined system which is expressed as 

/n^(r) ,Q „ fn'^(r) ,o „ f n^ (f)n'^(f') .,,•,, 1 f n^ (f)n^ (f') ,, ,o , 
^— (fr + Z / ^—^ (fr + E^- \' _^\ ' d^rd^r' + - ' ^ ' ^ ' -i^-^^^-' 
r . r . \r -r' \ 2 



The first two terms, Tg and Tp denote respectively, the kinetic energy of electrons and the positron, the next two terms 
represent the attractive and the repulsive interaction energies of the electrons and a positron with the nuclear charge. 
Ex is the exchange energy while the last two terms signify the electron-electron and electron-positron interaction 
energies respectively. The expressions for the calculation of exchange energy and electron-electron interaction energies 
can be found in the appendix of Ref. [g8| . The electron-positron energy can be obtained by following the steps exactly 
similar to the electron-electron case and emerges as 

J" 



E^-Y.11 ^"'™ /^' dr r'2 dr' RUt) R+^Hr')^ 

7 111 I " '~>. 

(II" I \ (I I" 

X (2Z + 1) (2r + l) 



TO -m/VO 

, I' I" V \ fV I" l'\ 
X , . , L „ J- (9) 



-to'/ VO 

Here, Rni and i?„/;/ are respectively the radial parts of the electron orbitals and the positron orbitals, and Nnim is 
the orbital occupancy. The 3j symbols in Eq. (^) arise due to integration over the solid angle Vl = ^{9, (j)). 

The Herman-Skillman code [^, modified for the Harbola-Sahni potential has been further modified in order to 
incorporate the positron. The calculations are carried out in the ceniraZ-field approximation for the systems Li~ 
through F" except for Be^ and N^ and the halide ions. In order to obtain the ground state of an anion A^, we start 
with converged potential of the neutral atom A and perform self-consistent calculation. The converged potential of 
the anionic system A" was then taken as a starting potential for the anion-positron self-consistent calculation. This 
was done in order to achieve fast convergence. In the following section we present our main results. 

III. RESULTS 

The total energies of the anion-positron bound states calculated in the present formalism and the corresponding 
restricted HF (RHF) energies for the positron in different states are displayed in Table |. The RHP numbers for total 
energies, positron and positronium affinities against which we compare our results are due to Patrick and Cade |l2] 
and Cade and Parazdel ||I^. The present total energies are in good agreement with the HP energies. The differences 
in parts per million between the energies of the present work and those of the HP theory are given in Table O for the 
anion-positron bound state (the positron is in the Is orbital). The differences diminish with the size of the anion. 
It is also evident from the table that the calculated total energies are slightly higher than the HP energies. This is 
expected since the HS orbitals differ from the HP orbitals which variationally minimize the total energy. 



The positron affinity is defined as 

P.A. = E{A-) + E{e+) - E{[A--e+]), (10) 

where E{[A^\ e+]) denotes the energy of the anion-positron bound system. Positive value of the P.A. indicates that 
iiJ[A~;e"'"] is a bound state, that is, the system A~ will bind a positron. The calculated positron affinities and 



the negative of the positron eigenvalue are given Table in [11 , also given are the corresponding HF values given for 



comparison. For all the systems, the PA is positive indicating their stability with respect to dissociation into an 
anion and a positron. The present values of the positron energy eigenvalues are in general higher in magnitude than 
the corresponding HF values and lie between the PA calculated from Eq. (HG) and the recent diffusion quantum 
Monte Carlo (QMC) values [|5[. The available diffusion Monte Carlo values [|5| for the positron affinities for different 
systems (in ground state) are 6.507(Li-), 6.015(5-), 5.941(C-), 5.862(0-) and 6.170(F-) in eV . 

Another binding energy of importance is the positronium affinity, a positive value of which means the system 
[A^ ] e^] is stable with respect to break up into the positronium {Ps) and a neutral atom A. The binding energies or 
positronium affinities can be computed in various ways Iq] . We compute the positronium affinities using the following 
two definitions: 

PsA = E{A)^^E{Ps)~E{[A-]e"^]) (11) 

and 

PsA = E.A. + P.A. + E{Ps). (12) 



Positronium affinities computed from Eq. (|l l[) are compared against the restricted Hartree-Fock (RHF) positronium 



affinities |l4,M in Tabic III . While for all the systems investigated the present and the RHF values of positronium 
affinities (calculated using Eq. (O)) arc in good agreement, no system is stable with respect to dissociation into a 
neutral atom and positronium: the positronium affinities for all systems are negative. 

In order to calculate the positronium affinity using Eq. ([L2|), we choose —e,nax of anionic system for the electron 
affinity (EA) as it is empirically found [|T| that in the exchange-only work formalism and the HF theory, the —tmax 
of the anionic system is, in general, a better estimate of EA than those obtained from the difference of self-consistent 
total energies of the atom and the corresponding anion. Further, it is observed that such estimates of EA within the 
present formalism are closer to the experimental EA [p2| than those obtained in HF theory by means of Koopmans' 



theorem ||30| . For positron affinity wc employ — e+ (third column in Table III) since this quantity is, in general, in 
better agreement with the accurate QMC positron affinity than the one obtained by taking difference of self-consistent 
energies (using Eq. ([lO|)). The positronium affinities thus calculated are also given in the last column of the Table 
m. These values of PsA are less negative than the PsA computed as the difference of the self-consistent energies 
(Eq.([ll|) leading in some cases to the binding. The systems [0~\e'^2, [F-;e+j,], [-F-;e^], [Cl~\e'^2 and [Br^;e'^^ 
are found to be stable against the dissociation into the positronium and an atom. This binding may be attributed to 
the accurate asymptotic structure of the work formalism HF approximation. 



We finally present the one electron properties such as < r" > expectation values for halide anions in the Table |rV . 
The one electron expectation values are in good agreement with their HF counterparts. The computed < l/r> values 
are slightly larger than the HF values, implying the slight increase in the positron density towards the nucleus which 
therefore, should be compensated by small reduction in the long-range of the positron density leading to smaller < r > 



and < r^ > . This is indeed the case as can be seen from the Table IV . 

We have also investigated the binding of positron to neutral atoms in the spirit of Patrick and Cade, by starting 
with the anion-positron bound state and reducing the ionicity of the system to obtain desired neutral system. It was 
found that the neutral-atom positron bound state does not exist in the exchange-only work formalism. This, however, 
is not surprising as the present treatment lacks the electron-positron and electron-electron correlation effects which 



are crucial to permit such a binding |gl[. Further, it has been rigorously shown by Pathak JS^] that the deviation 
from the spherical symmetry is a necessary condition in order that the positron binding to neutral atom would occur. 
The present treatment can be extended to include the correlations effects, namely the electron-electron Coulomb 
correlations, correlation contribution to the kinetic energy and the electron-positron correlation. The first one can be 
incorporated by modeling the correlation second-order density matrix as suggested by Levy and Perdew [g^ while 
the second one can be derived in terms of density matrices via virial theorem |25| . The electron-positron correlation 
potential l^^l can be added in an ad hoc manner to the effective potential. 

IV. CONCLUSIONS 

In the present work, the positron binding to negative ions is investigated within the exchange-only work formalism. 
The work formalism of Harbola-Sahni seems to provide the Hartree-Fock level description of the electron-positron 
system as can be seen from the agreement between the present values of positron and positronium affinities and 
their restricted Hartree-Fock counterpart. The advantage of the work formalism is that its effective potential is local, 
orbital independent and therefore computationally cheaper. The systems [0~;e^^], [F~;e^g], [F~]e2p], [Cl~;e'^^] and 
[Br^ ; ejj.] are found to be stable against the dissociation into the positronium and the corresponding atom. 



TABLE I. Negative total energies (in Hartree a.u.) of [A ; e^] bound system calculated in central field approximation 
within the work formalism. Also given are the total energies in HF theory. 



System 



Work formalism 



HF 



System 



Work formalism 
-E 



HF 

-E 



Li" Is^ 2s^ 


7.4270 


7.4282 


Li^ : Is 


7.5286 


7.5299 


Li' : 2s 


7.4748 


7.4760 


Li- : 3s 


7.4528 


- 


Li- : 2p 


7.5017 


7.5030 


Li- : 3p 


7.4653 


- 


Li- : 3d 


7.4752 


7.4765 


Li Is^ 2s^ 


7.4316 


7.4328 


B- ls^2s^2p' 


24.5156 


24.5192 


B- : Is 


24.6495 


24.6531 


B- : 2s 


24.5733 


24.5769 


B- : 2p 


24.6108 


24.6202 


B- :3s 


24.5477 


- 


B- :3p 


24.5610 


- 


B- :3d 


24.5694 


24.5757 


B ls^2s^2p^ 


24.5261 


24.5292 


C- ls'^2s^2p^ 


37.7041 


37.7088 


C- : Is 


37.8563 


37.8610 


C- : 2s 


37.7671 


37.7718 


C- : 2p 


37.8040 


37.8087 


C- : 3s 


37.7384 


- 


C- : 3p 


37.7513 


- 


C- :3d 


37.7584 


37.7632 


C ls^2s^2p^ 


37.6847 


37.6887 


Q- ls^2s^2/ 


74.7849 


74.7897 


Q- : Is 


74.9583 


74.9630 


Q- : 2s 


74.8534 


74.8582 


Q- :2p 


74.8940 


74.9026 


Q- :3s 


74.8214 


- 


Q- :3p 


74.8350 


- 


Q- :3d 


74.8403 


74.8461 


O ls^2s^2p'' 


74.8050 


74.8095 



F- Is22s22/ 


99.4543 


99.4594 


F- : Is 


99.6383 


99.6434 


F- :2s 


99.5253 


99.5305 


F- :2p 


99.5641 


99.5692 


F- :3s 


99.4917 


- 


F- :3p 


99.5048 


- 


F- :3d 


99.5095 


99.5147 


F ls^2s^2p^ 


99.4046 


99.4095 


cr [Ne]3s^3p'^ 


459.5640 


459.5769 


cr : Is 


459.7071 


459.7189 


cr : 2s 


459.6243 


459.6373 


cr : 2p 


459.6625 


459.6754 


cr : 3s 


459.5972 


- 


cr : 3p 


459.6107 


- 


CI [Ne]3s''3p^ 


459.4697 


459.4830 


Br- [Ar]4s^3d^°4p'^ 


2572.523 


2572.5363 


Br- : Is 


2572.656 


2572.6695 


Br- : 2s 


2572.5803 


- 


Br- : 2p 


2572.6177 


2572.6311 


Br [Ar]4s^3d^°4p'^ 


2572.229 


- 



TABLE II. The total ground state energy differences between tfie present work and the HF theory in parts per million for 
anion-positron bound states (the positron is in the Is orbital). 

system differences(ppm) 

[Li-,e+] 173 

[B-,e+] 146 

[C-,eU 124 

[0-,e+] 59 

[F-,e+] 51 

[Cr,e+] 26 

[Br-,e+] 5 



TABLE III. Positron eigenvalues, affinities and positronium affinities in eV 



system 


n+l+ 


— f 


; + 


PA 


PsA 


PsAt 






HS 


RHF 


HS 


RHF 


HS 


RHF 




Li- 


: Is 


3.049 


2.996 


2.765 


2.766 


-4.161 


-4.159 


-3.341 


Li- 


: 2s 


1.347 


1.329 


1.306 


1.301 


-5.625 


-5.624 


-5.0427 


Li- 


: 3s 


0.770 


- 


0.770 


- 


-6.169 


- 


-5.620 


Li- 


:3p 


1.064 


- 


1.042 


- 


-5.883 


- 


-5.326 


Li- 


: 3d 


1.339 


1.331 


1.312 


1.312 


-5.614 


-5.613 


-5.051 


B- 


: Is 


3.785 


3.778 


3.644 


3.642 


-3.422 


-3.428 


-2.395 


B- 


: 2s 


1.584 


1.582 


1.714 


1.569 


-5.372 


-5.502 


-4.596 


B- 


:2p 


2.557 


2.799 


2.590 


2.748 


-4.240 


-4.323 


-3.623 


B- 


: 3s 


0.876 


- 


0.875 


- 


- 


- 


-5.304 


B- 


:3p 


1.212 


- 


1.235 


- 


- 


- 


-4.968 


B- 


: 3d 


1.442 


1.540 


1.464 


1.536 


-5.450 


-5.535 


-4.738 


c- 


: Is 


4.204 


4.218 


4.142 


4.141 


-2.131 


-2.21 


-0.854 


c- 


: 2s 


1.712 


1.718 


1.714 


1.713 


-4.558 


-4.540 


-3.346 


c- 


:2p 


2.712 


2.135 


2.718 


2.718 


-3.554 


-3.535 


-2.346 


c- 


: 3s 


0.931 


- 


0.933 


- 


-5.339 


- 


-4.127 


c- 


:3p 


1.282 


- 


1.284 


- 


-4.988 


- 


-3.776 


c- 


: 3d 


1.478 


1.480 


1.478 


1.480 


-4.795 


-4.773 


-3.580 


o- 


: Is 


4.784 


4.769 


4.718 


4.716 


-2.628 


-2.662 


1.443 


o- 


: 2s 


1.869 


1.865 


1.864 


1.862 


-5.483 


-5.475 


-1.472 


o- 


:2p 


2.917 


3.079 


2.969 


3.070 


-4.378 


-4.268 


-0.424 


o- 


: 3s 


0.993 


- 


0.993 


- 


-6.354 


- 


-2.348 


o- 


:3p 


1.350 


- 


1.363 


- 


-5.984 


- 


-1.991 


o- 


: 3d 


1.498 


1.532 


1.508 


1.532 


-5.839 


-5.805 


-1.843 


F- 


: Is 


5.061 


5.048 


5.007 


5.006 


-0.441 


-0.434 


3.126 


F- 


: 2s 


1.932 


1.936 


1.932 


1.933 


-3.516 


-3.506 


0.003 


F- 


:2p 


2.993 


2.992 


2.998 


2.987 


-2.460 


-2.452 


1.058 


F- 


:3s 


1.020 


- 


1.018 


- 


-4.430 


- 


-0.015 


F- 


:3p 


1.377 


- 


1.374 


- 


-4.073 


- 


-0.598 


F- 


:3d 


1.503 


1.502 


1.503 


1.503 


-3.946 


-3.937 


-0.432 


cr 


: Is 


3.928 


3.922 


3.948 


3.894 


-0.340 


-0.350 


1.052 


cr 


: 2s 


1.645 


1.644 


1.641 


1.642 


-2.590 


-2.600 


-1.231 


cr 


:2p 


2.687 


2.687 


2.680 


2.680 


-1.554 


-1.560 


-0.189 


Br- 


: Is 


3.655 


3.653 


3.619 


3.626 


-0.620 


-0.690 


0.409 


Br- 


:2s 


1.568 


- 


1.559 


- 


-2.680 


- 


-1.678 


Br- 


:2p 


2.585 


2.588 


2.577 


2.577 


-1.670 


-1.740 


-0.661 



t Calculated using Eq. ([12) (see text for details) 



TABLE IV. The radial expectation values < i >, < r > and < r^ > for the Is positron orbital in halide ions calculated 
within the work formalism and the corresponding Hartree-Fock values (the values are in Hartree atomic units). 



system 




< 


i> 






<r> 






<r' 


> 






HS 






HF 


HS 




HF 


HS 






HF 


[F-;e+] 


0.2959 






0.2948 


4.067 




4.080 


19.454 






19.572 


[Cr;e+] 


0.2194 






0.2189 


5.2987 




5.3085 


32.2031 






32.3131 


[Br-;e+] 


0.2018 






0.2015 


5.7178 




5.7198 


37.1800 






37.2715 
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